Comparison Results for Garch Processes 

Fabio Bellini * Franco Pellerey } Carlo Sgarra^and Salimeh Yasaei Sekeh 

March 2, 2013 



Abstract 

We consider the problem of stochastic comparison of general Garch- 
like processes, for different parameters and different distributions of the 
innovations. We identify several stochastic orders that are propagated 
from the innovations to the Garch process itself, and discuss their inter- 
pretations. We focus on the convex order and show that in the case of 
symmetric innovations it is also propagated to the cumulated sums of the 
Garch process. More generally, we discuss multivariate comparison results 
related to the multivariate convex and supermodular order. Finally we 
discuss ordering with respect to the parameters in the Garch (1,1) case. 
Keywords: Garch, Convex Order, Peakedness, Kurtosis, Supermodularity. 



1 Introduction 

An extensive literature is available on applications of stochastic orders to finance 
and insurance markets. The implications of stochastic orders for derivative pric- 
ing and risk management are relevant. The increasing dependence of european 
option prices by the riskiness of the underlying it is a well known property for 
basic models like that of Black- Mcrton-Scholes, in which riskiness is expressed in 
terms of the logreturns distribution variance: the uncertainty is quantified there 
through the dispersion around the expected value and the distribution functions 
can be ordered according to their "peakedness"; the larger is the dispersions, 
the higher the option prices. This very elementary and intuitive observation for 
simple models become more involved when turning attention to more complex 
models, where a more rigorous approach is necessary in order to avoid wrong 
conclusions. 

The Black-Merton-Scholes model is nowadays considered fairly inadequate to 
describe the asset price dynamics; several empirical facts cannot be explained on 
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the basis of this model: some statistical features exhibited by logreturns like fat 
tails, volatility clustering, aggregational Gaussianity and the so-called leverage 
effect are completely outside of the prevision properties of the Black-Merton- 
Scholes model. Moreover a very relevant phenomenon exhibited by option prices, 
the "volatility smile" (and its term structure) cannot be explained on this model 
basis. In order to provide a more satisfactory description several different models 
have been introduced. Some of these models introduce a stocastic dependence 
in volatility and/or jumps in asset logreturns (and/or in volatility) dynamics 
both in continuous and discrete time setting. 

Among discrete time models introduced in order to remove some of the 
Black-Merton-Scholes model drawbacks, the class of Autoregressive Conditioned 
Heteroschedastic (ARCH) models introduced by Engle [6] and their general 
extension (Garch models) proposed by Bollerslev in [4] have risen considerable 
interest. 

Several results related to stochastic orders are available for the continuous 
time models class: in [2], where a systematic investigation on semimartingale 
models is performed; the models considered there include the Heston and the 
Barndorff-Nielsen and Shephard models. In [IT] T. M0ller provides some re- 
sults on stochastic orders in a dynamic reinsurance market where the traded 
risk process is driven by a compound Poisson process and the claim amount 
is unbounded. Stochastic order properties have been used to obtain bounds 
for option prices in incomplete markets; the literature focused on this subject 
is quite extensive and we just mention the papers by El Karoui et al. [5], by 
Bellamy and Jeanblanc, and by Gushchin and Mordecki [8]. 

The purpose of the present paper is to present a systematic investigation of 
stochastic orders propagation in a Garch context. 

Comparison with stochastic orders in incomplete market models can give 
rise to different classes of problems: first can be considered the comparison of 
models under the same probability measure but with different parametric spec- 
ification, second it can be examined the problem of comparing the same model 
under different probability measures; as a matter of fact, when markets are 
incomplete, there are several probability measures equivalent to the historical 
one, under which the dynamics of prices can be given. In this paper we shall 
focus on the first class of problems mentioned: we shall provide a systematic 
comparison of logreturns and then of prices when the model parameters change, 
but the dynamics is specified under the same probability measure. In a Garch 
context the parameters entering into play are three parameters assuming a nu- 
merical value and the innovations, which are random IID variables for which 
the density function is assigned. We just mentioned that stochastic order re- 
sults have important implications on option pricing and this holds true also in 
a Garch context: in particular, convex order relations on logreturn sums imply 
increasing convex order relations on the underlying price, hence on european 
call option prices and this can be considered the main relevance of our results 
from an application viewpoint. 

We like to present a numerical illustration in order to motivate our investiga- 
tion. Fig 1 compare the densities of the logreturn sums in a Garch (1,1) model 
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Figure 1: Comparison of Logrcturn Sum Densities 



with respect to variations in the parameters ao,ai,/3i; in the continuous line 
all parameters assume the value 0.2, while the other lines represent the same 
density, but with parameters a = 0.5,ai = 0.5 and /?i = 0.5 respectively (the 
parameter assuming the value 0.5 in indicated by a capital letter in each curve 
caption). The sums include the first 50 terms of the logreturn sequence. 

The innovations, assumed to be independent and identically distributed, 
and the initial variance, are chosen as standard Gaussian random variables. It 
is evident how the usual real numbers order relation between parameters implies 
some ordering on logreturn sums; it then arises as a natural question to ask if 
this simple remark can be made rigorous and if this conclusion can be cast into 
a more general framework considering different kind of stochastic orders and 
comparing stochastic order propagation from innovations to logreturns and to 
logreturns sums. 

In the following section we present briefly the Garch models and an auxiliary 
lemma. In Section 3 we present the univariate stochastic comparisons for logre- 
turns in a Garch setting. In Section 4 we provide results on some related orders, 
while in Section 5 the convex order propagation of the logreturn sums is inves- 
tigated together with its implications on asset and european call option prices, 
while their multivariate convex order propagation is the subject of Section 6 . 
We focus our attention on the Garch(l,l) case in the last section. 

The main results provided in this paper require the assumption of a sym- 
metric probability density for innovations. The extension of the stochastic order 
propagation analysis presented here to the case of non-symmetric innovation 
densities would be of great interest: both the convex order propagation result 
and the comparison with respect to parameter variations will be the subject of 
future investigation in this more general setting. Moreover the identification 
of some convex multivariate order which naturally propagate from innovations 



to logreturns is another target on which our research interest will be focused. 
These will be the subject of our future work. 



2 General Garch models 

We consider Garch models of two different very general forms; the first model 
(Ml) is: 

X n = (T n s n , n = 0, 1, .. 
e n _L a n , E[e n ] = (1) 

Cn+l = f{\£n\ ,CT„) 

with f 1 : — > increasing and componentwise convex (ccx for brevity). 
The second model (M2) is 

X n = cr n e n , n = 0, 1, .. 
e„ J- a n , E[e n ] = (2) 

with f 11 : — >• M + increasing and ccx. In both cases the innovations e„ are 
independent and identically distributed (i.i.d.). When discussing the propaga- 
tion of variability orderings, the normalization requirement -E[e^] = 1 will be 
dropped. The difference between model Ml and model M2 is that in the first 
case the recursive dynamics is defined in terms of the volatility a n , while in the 
second it is defined in terms of the variance cr^ . 

The usual Garch (1,1) model is a particular case of both Ml and M2, and is 
defined as follows: 

X n = a n e n , n = 0, 1 
n+i =a + a 1 X? l + (3 1 a, 

with ao,ai,Pi > and a\ + fi\ < 1, in order to guarantee covariance station- 
arity. Both models start with a possibly random o~o > 0, by drawing a random 

£ . 

The recursive equations for OYi+i and er^ +1 are examples of "stochastic recur- 
rences" in the sense of Chapter 4 of [TU]. For the explicit expression of their 
solutions, we introduce the following notations: 

<r n +i ■= gi+i(vo, M , \s n \) (4) 

2 J I / 2 Jl _2\ 

°n+l ■— i/n+lV cr 0i fc 0i •") t n) 

As in [10] , we have the following: 

Lemma 1 Let g^ +1 , g^+i '■ R™ +2 — > R+ be defined as in Then g^ +1 and 
gn+i are increasing and componentwise convex. 

Proof. We have clearly g£ +1 (<7 , Nl > ••-> l e n|) = /'(M , / J (ko| , 0o)) 
and 1 (a 2 , si , £ 2 n ) = f n (e 2 n , . .. , f n (eg , a 2 ) ) ; since / J and / IJ arc increasing 
and ccx, then also their compositions in Q are increasing and ccx (see Meester 
and Shanthikumar [9] for this last assertion and further properties of increasing 
and ccx functions). ■ 



2 , fl _2 (3) 
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3 Univariate comparisons of X, 



The aim of this section is to establish comparison results for X n when the 
distributions of the innovations are changed from Ek to Ek ■ In order to establish 
these results, the assumption that the innovations are identically distributed is 
not necessary (while the independency assumption is essential) . In the following 
theorems only the distribution of a single innovation e k will be changed, and 
the impact of this change on X n will be investigated. 

We recall the definitions of the basic stochastic orderings: 

Definition 2 The random variable X is said to be smaller than Y in the usual 
stochastic order [convex order, increasing convex order], denoted by X < st Y 
[X < cx Y, X <i CX Y], if E[(fi(X)] < E[cf)(Y)] for all increasing [convex, increas- 
ing convex] functions <fi '■ K — > K for which the expectations exist. 

We will see that in the general context of models Ml and M2 the orderings 
that are naturally propagated from the innovations £k to X n are the < s t and 
the <i cx ordering between absolute values or squared variables. This clearly 
completely modifies their interpretation; in particular, in the next section we will 
see that the < st ordering between absolute values or squares can be interpreted 
as a variability ordering, while the <i cx ordering between absolute values or 
squares can be interpreted as a kurtosis ordering. 

In order to establish these results, we proceed in two steps: first we consider 
the volatilities o~ n and then the variables X n . The first step is an immediate 
consequence of Lemma [TJ 

Theorem 3 Comparisons of a n and a\ 

a) Let a n+ i be as in ^j) and \eu\ zf> s t |£fc|,' it follows that cr n +i ^ s t &n+i- 

b) Let a n+ i be as in (QP and \e k \ < icx |e fc |; it follows that a n+1 < icx <J n+ i. 

c) Let a^+i be as in (Hp and e\ ~< st ej,; it follows that c^+i < s t o^+i- 

d) Let <t^ +1 be as in and e\ ~<i CX ef ; it follows that -<i CX &n+i- 

Proof. Since from Lemma [T] in model Ml <r„+i = <7^ +1 (<to, |eo| , |^n.|) 
with <?„_|_ x increasing and ccx, item a) and b) follows respectively fromTheorem 
1.A.3 in [T2] and Theorem 4. A. 15 in [T3]. Similarly, since from Lemma [T] in 
model M2 cr, 2 l+1 = 5^+i(cq, e§, e 2 n ) with g^+i increasing and ccx, from the 
same theorems we get c) and d). ■ 

The comparison results for a n and <r 2 lead to the following comparisons of 
the variables X n : 

Theorem 4 Comparisons of X n 

a) Let X n be as in (Qp and \ek\ ^ s t |ejfe|; it follows that \X n \ ^ st 

b) Let X n be as in {TJj and \sk\ diicx it follows that \X n \ <i C 

c) Let X n be as in (0) and e\ -< st it follows that X% < st X 2 . 

d) Let X n be as in |2[) and si <i CX e|; it follows that X 2 <i CX X 2 
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Proof. Since \X n \ — <j n \s n | and X 2 = a^e^, w ith °n independent from e n , 
items a), and c) follow from Theorem 1.A.3 (b) in |13) . Similarly item b) and 
d) follow from Theorem 4. A. 15 in [13]. ■ 

A natural question that arises at this point is if also the convex order is 
propagated, that is if Ek d:cx £fc X n -< cx X n . We prove that this is indeed 
the case for model Ml. We start with a simple lemma: 

Lemma 5 Let a and a be nonnegative, with a -< st a. Let e be independent 
from a and a, with E [e] = 0; then as -< cx as. 

Proof. We can construct identically distributed copies of a and a on the 
same probability space, such that a < a a.s. Then for each realization of er, we 
have that ere < cx ae. From Theorem 3. A. 12 in [13] it follows that ae < cx ae. ■ 

Theorem 6 Propagation of convex order 

Let X n be as in (QJ) and < cx it follows that X n ^< cx X n 

Proof. First of all we remark that since £k dicx £k, it follows that \ek\ diicx 
|efc|- Indeed, for each <f> increasing and convex, the composition 0(|...|) is convex; 
this implies that 23[0(|efc|)] < i?[0(|efc|)], that is \Ek\ diicx |efc|- From Proposition 
[3jitem b), it then follows that a n +i diicx 5v»+i- From Theorem 4. A. 6 in [13] there 
exists a random variable ~a n +\ with <7 n +i ;K s4 ~5 n +i ^cx &n+i- By Lemma [5] 
we have that c^+i -< s t ~a n +i implies that a n+ ie n +i ^ cx a n +i£ n +i] on the other 
hand a n +i dicx a n +i implies that a n+1 e n+ i < cx a n+ i£ n+ i. By transitivity we 

get Cn+iEn+i < cx <7 n +i£ n +i. ■ 

4 The relevant orderings 

In the preceding section the orderings defined by \X\ ^< st \Y\, X 2 ^< st Y 2 , 
\X\ diicx \Y\, X 2 <i CX Y 2 have arisen naturally. 

In order to better understand their meaning, in the following lemmas we identify 
some necessary and sufficient conditions in the continuous and symmetric case. 
We have the following: 

Lemma 7 Let X and Y be symmetric with continuous distributions F and G. 
The following conditions are equivalent: 

a) X 2 < st Y 2 ; 

b) \x\ ± st m- 

c) X ^peak Y , where -< pea k is the peakedness ordering introduced in 

d) S~ (G — F) = 1 with sign sequence +, — , where S~ (G — F) is the number of 
intersections between G and F as defined in (1.A.18) of 

Proof. The equivalence of a) and b) is an immediate consequence of Theo- 
rem 1.A.3 in [T3]. The equivalence of b) and c) is the definition of the peakedness 
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ordering, while the equivalence between c) and d) follows from Theorem 3.D.1 
in [13]. ■ 

Lemma 8 Let X and Y be symmetric with continuous distributions F and G. 
The following conditions are equivalent: 

a) X 2 ^ x Y 2 

b) 00 F( u ) ■ u ■ du < °° G(u) ■ u ■ du for each x > 0, where F(u) = 1 — F(u) 
and G{u) = 1 - G{u) 

c) E[(X 2 - k)+] < E[(Y 2 - k)+] for each k > 

Proof. Under our hypothesis F x *(t) = 2F(v / *)-l andF x2 (i) = 2-2F(y/t), 
for t > 0. The equivalence of a) and b) follows then from Theorem 4. A. 2 in [13] 
with a simple change of variable. The equivalence between a) and c) is also a 
consequence of Theorem 1.5.7 in [TU] ■ 

The first lemma shows that for symmetric variables the orderings \X\ ^< st 
\Y\ and X 2 ^< st Y 2 are variability comparisons equivalent to the peakedness 
ordering, that in this case boils down to item d), that is the validity of a single 
cut condition between the distribution functions. In the typical econometric 
applications these orderings are however not very relevant since the innovations 
satisfy E[e\] — 1, and hence E[e 2 k ] ^ st E[ef,] would imply e| — st e|. 

In the normalized case the ordering X 2 <i CX Y 2 becomes equivalent to 
X 2 -< cx Y 2 ; we prove a sufficient and a necessary condition for it. 

Lemma 9 Let X and Y be symmetric with continuous distributions F and G 
and with E[X 2 ] = E[Y 2 } = 1. 

a) If the densities of X and Y cross 4 times, with the density of X being lower in 
the tails and in the center, and higher in the intermediate region, then X 2 ^i CX 
Y 2 . 

b) If X 2 ~<i cx Y 2 and X andY have finite fourth moments, then ^(X) < faiY), 
where $2 is Pearson's kurtosis coefficient. 

Proof, a) Under our hypothesis f x *{t) = for t > 0. Since X and Y 

are symmetrical, we have that the four intersection points between the densities 
/ and g are symmetrical with respect to the origin. Hence the densities of X 2 
and Y 2 cross in two points and since E[X 2 } = E[Y 2 } from Theorem 3. A. 44 in 
[T3] we have that X 2 < cx Y 2 . 

b) In our case ^(X) = E[X 4 } and hence the thesis follows from the definition 
of the convex order. ■ 

This lemma shows that the comparison X 2 -<i CX Y 2 can be interpreted 
as a classical kurtosis ordering; the cut condition is usually referred in the 
kurtosis ordering literature as a Dyson-Finucan condition (see for example [7 
for a review). 
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5 Convex comparisons for total logreturns 



In financial applications the variables X n typically represent logreturns, that are 

n 

additive quantities. The over-the-period total return is given by S n — X n - 

k=i 

It is therefore natural to ask if some of the comparison results of Section 2 do 
extend to the variables S n . In this section we consider the case of the convex 
order, that is, we wonder if eu dicx £k => S n < cx S n . The problem is not trivial 
since S n cannot be expressed as a sum of independent variables, so the standard 
results about convex ordering of sums cannot be applied; we are able to prove 
a positive result in the case of model Ml and for symmetric innovations. We 
start with a basic lemma: 

Lemma 10 Let <p € C 2 (K) be convex and gt G C 2 (M) be convex and nonnega- 
tive. Let a, b G R and P m := { — 1, l} m . It follows that 

h(u) = \^ (f> I a + bu + y^Pigi(u) 



/ J Y 1 1 1 / , 1 

peP m \ i=i 



is convex. 

Proof. We can compute 



h'(u)= <t>' {a + bu + J2 Pl g t (u)\ ■ ib + ^p^u) 

peP m \ *=i / \ »=i j 



h"(u) = £ 4>" (a + bu + ^pMu) ■ U + Y J Pi9» + 

pGPm V i=l / V i=l / 

(m \ / m \ 

a + bu + ^p l g t {u) J • I ^Pig" (u) J 
»=i / \*=i / 

The first term is positive; the second is given by 

(m \ / m \ 

a + bu + ^ Pi9i («) J ' ( X] ^-9" ( u ) J 
i=l / \*=1 / 

Let us denote with P a random vector with a discrete uniform distribution on 
P TO ; clearly E\P] = 0, the components of P are independent and 

A m = 2 m E [</>' (a + bu + g(u) ■ P) (g» ■ P)] 

Since the functions <p' (a + bu + g(u) ■ P) and g"(u) ■ P are componentwise in- 
creasing, from the covariance inequality it follows that 

A m = 2 m E[^{g{u)-P){g['(u)-P)]> 

> 2 m E[<p'(g(u)-P)]E[(g;'(u)-P)]=0 
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We remark that in this lemma the smoothness requirements on <f> and on the 
gi can be dropped; we preferred this formulation in order to simplify the proof. 
Since in this section we consider only model Ml, we define 

g n (a ,£ , ...,£„_i) := g n (o- , |e | , |£n-i|); 
from Lemma [1] it is clear that g n is even and ccx. We have 

S n = Xq + Xi + ... + X n = a e + o\E\ + ... + a n e n = 

CTO^O + 5l(co,£o)£l + ■•■ + 0n(o'O>£Oj -,£n-l)£n = (5) 
Sni^O, £ 0: •■•7 Sn) 

The main problem in proving the propagation of convexity to the sums is that S n 
is not a ccx functions of the innovations Sk', indeed, each in (J5J) is multiplied 
by a possibly negative innovation £& . This prevents the applications of standard 
results and requires the development of a specific technique based on Lemma 
[TU1 The basic idea is that in the case of symmetric innovations it is possible 
to restore the convexity by averaging over all the possible sign changes, as in 
Lemma [TU1 This will be done in a recursive way; we start with the following: 

Lemma 11 Let X n and S n be as in (QP and Let <p be convex and Ei be 

symmetric. Then the function 

h(a ,e Q , e k ) := E ek+1> ... >Sn [4>(S n (a Q ,SQ, £„))] (6) 

is convex in Ek for each fixed value of ao, sq, £k-i- 

Proof. To avoid notational burdening we drop the arguments of the func- 
tions gi. Since the innovations are symmetric and gi is even, we can write 

E ek+1 ,..., £n [(f>(Sr l (ao,eo, ... ,£„))] = E eh+u ... t e n [<t>(a eo + ... +#„£„)] = 

E Sk+ll ...,e n [ ^ ^("OEO + ■■■ +Pl3fc+l£fc+l + ... +p„-fc.9„£„)l{ £fc + 1 >o,..., £ „>0}]- 

Denoting by 

h(<ra,eo, ... ,£fc, — ,e n ) = 2J 4>(< J £ 0+9l£l+---+9k£k+Pigk+l£k+l+---+Pn-kgn£n), 

we have that 

h{o~o,£o, £fc-l, £fc) = -Ee fc + 1 ,...,e Il [l{e fc + 1 >0,...,e„>0}^( <J 0,£07 •■•,£fc, •■•,£«)] 

and h is convex in Sk from Lemma [TOl It follows that also h(o~o, £o, £fe-i, £fc) 
is convex in eu for each value of ao, £q, £fe_i. ■ 

We can finally state the result on the propagation of the convex order to S n : 
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Theorem 12 Let X n and S n be as in (QJ) and 0). Let Ej be symmetric. If 
also e k is symmetric and e k > cx e k , then S n := S n (ao,£o, .,e k , ...,£„) > C x 
(cr ,eo, -,£fc, •■•,£«)• 

Proof. Let (?!> be convex. From the independence of the Ei we can write 

E[4>(S n )} = E £0t ,,, t£k _ 1 E £k E £k+lt ,,, t£n [(j)(S n (a ,e , ...,e k -i,e k ,e k+ i, ...,e„))] = 
= E £0t ,,, t£k l E £k [h{o- f)l eo, ...,ek-i,£k)] 

where as in ([6]) 

h(<Jo, £o, £fc-l, £fc) ■— ^e fc + 1 ,...,e„ [4>(S n ((J Q , £ , £fc-l, £fc, £fe+l, ...,£„))] 

is a convex function of e k for each value of o~$, £o, £fc-i from Lemma [TT1 
It follows that 

Ee k [h( a 0> e 0, £fc-i, £fc)] > E £k [h(a , £o, £fc-i> £fc)] 

that gives 

E[cj)(S n )} = E £(u ..., £k _ 1 E £k [h(o-o,eo,...,Ek-i,e k )}> 

> ^e o ,...,e fc -i- E e fc [ /l ( cr 0,£0, ■•■,£fc-l, £fc)] = ^(^n)] 

that is S n ^_ C x S n . ■ 

6 Multivariate comparisons of logreturns 

Until now we have been considering only univariate orderings of the X n . Having 
established also a convex comparison result for the sums S n , it is natural to 
wonder whether more general multivariate comparisons do hold, that is if we 
can prove that e k <cx £k => (Xi, X n ) ■< (Xi, X n ) for some multivariate 
convexity ordering ^ to be precise in the sequel. Before stating a positive result, 
we recall two basic definitions: 

Definition 13 A function ip : 1™ — > M. is directionally convex if for any 
Xi, ...,X4 6 M™, such that xi < X2,X3 < X4 and xi +X4 = x 2 +X3, it holds that 

</>(x 2 ) + ^( x 3) < ¥>(xi) + y(x4) 

Definition 14 A function ip : M™ — > M is supermodular if for any x, y G K™ 

satisfies: 

¥>(x) + ^(y) < ¥>( x A y) + (ys(x V y), 

where the operators A and V denote respectively coordinatewise minimum and 
maximum (see Section 7. A. 8 of JlSf ) 
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In the univariate case directionally convexity is equivalent to convexity, while 
in the multivariate case there are no implications between the two concepts. 
Directionally convexity implies supermodularity and its equivalent to super- 
modularity plus componentwise convexity. For smooth functions, directionally 
convexity is equivalent to the nonnegativity of all entries in the Hessian matrix, 
while supermodularity is equivalent to the the nonnegativity of all entries out of 
the principal diagonal. Clearly no implications exist between this concept and 
the usual convexity of </>, that corresponds to the positive semidefinitcness of 
the Hessian matrix. However, cf> is directionally convex and convex if and only 
if it is supermodular and convex. Finally, in the smooth case, tp is directionally 
convex if and only if its gradient is increasing. 

In order to establish multivariate comparison results, we need a generaliza- 
tion of Lemma [TOl 

Lemma 15 Let (f> G C 2 (K") be convex and supermodular and gi £ C 2 (R) be 
convex and nonnegative. Let P n := { — 1, 1}™. Lt follows that 

K u ) = < Mpi5i( u )'---'-P«5«( m )) 

is convex. 

Proof. Let's denote by yi the arguments of the function (f>; we can write: 

h'(u)= ^2 ^Pi9i(u)^-(pigi(u),...,p rn g m (u)) 

P eP m i=i ^ l 



h "( u ) = ^*^Pi9Ku)^—{pigi{u),.-.,p m g m {u)) + 

pep m i=i Vl 

+ ^2^2 PlP 3 9i( u )9j(u) g g (pigi(u),..., p m g m (u) )] 
The second term in the square brackets can be written as 







2, 



d 



2, 



dyidy. 



•(••) 



(Pg'(u)Y >0 



y^y^PiPjg'Mg'jiu) (pm(u), -,p m g m {u)) = Pg'(u) 

and is positive since the Hessian of <p is positive semidefmitc. 
The first term can be written as 

5DftSi'(«)^(Pifli(«). = 2 m E [V4> (pi5i(u), ..., P n9n(u)) ■ (iV'(tt))] 
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From the hypothesis on <f> all the components of V</> are increasing in pi, hence 
from the multivariate covariance inequality 

E [V<t>(pm(u),..., Pn g n (u)) ■ (Pg»)] > 

■ 

As in Lemma I10| the smoothness requirements on <p and gi can be dropped, 
but we added them in order to simplify the proof. The multivariate analogue 
of Lemma [TT1 is the following: 

Lemma 16 Let X n and S n be as in |7J) and Let <f> : R n+1 —> R be super- 
modular and convex and e, be symmetric. Then the function 

h k {x) := E[(j>(X ,...,X n )\e k = x] ; 



is convex. 

Proof. From the symmetry of the innovations we can write 

hk(x) = E £o ... £k _ l£k+1 £n [<j){<To£o,9i£i, -,gkX, , ...,g n e n )} = 

= E £o ... £k _ l£k+1 £n [l {£ >o} ^2 ^(o-qPoSo, giPi£i, gkPkX, g n p n £ n )} 

Since each gi is convex in e k , from Lemma 1151 it follows that for each o~o > 
and £i > 0, i ^ k, the function 

hk{x) = ^2 c/)(crop eo,g 1 p 1 e 1 ,...,g k pkX,...,g n p n e n ) 

pen, 

is convex. Averaging with respect to £j, with i ^ k, it follows that also h k {x) is 
convex. ■ 

This enables us to state our main multivariate comparison result: 

Theorem 17 Let X n and S n be as in (Qp and |5p. Let the Si be symmetric. If 
also e k *s symmetric and e k > C x £k, then, 

E [0(A O , X k , X n )] < E y>{X , X k , X n ) 

for every function (f> : R n+1 R supermodular and convex. 
Proof. From the previous lemma we have that 

E[<j>{X ,...,Xk,...,X n )] =E Sh [hk(e k )} <E £k [h k {I k )} = E^(X ,...X k ,...,X n ) 
■ 

We remark that we are not able to prove supermodularity or componentwise 
convex ordering of (Xq, ...,X n ); at the moment both hypothesys on <j> seem to 
be necessary for Lemma [T5l 
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7 The Garch (1,1) case 



We focus now on the Garch (1,1) model specified by 



'n+l 



X n 
■■ a 



^2 



aiX n + fcai 



(7) 



with ao,ai,Pi > and ol\ + j3\ < 1. For this model the recursive dynamic 
of the volatility or of the variance ((4]) can easily be explicitated as follows (see 



Tt + l 



zn^ 

fc=l i=l 



i + ai£„_ i+ i) 



(8) 



From this expression it is immediate that cr n+1 an d &n+i are nondecreasing 
functions of the parameters ao, aiand fi\ .We already remarked that this model 
is a special case of both Ml and M2, so all the comparison result for varying 
innovations of the preceding sections do hold. In this section we are interested 
in establishing comparison results for different parameters ao, ol\, j3±. As men- 
tioned in the Introduction, it is natural that an increase in ao, ot\, fi\ should 
correspond to an increase in the variability of X n and S n ; in this section we 
prove it rigorously. Without any additional effort, we can consider stochastic 
parameters ao, a±, f3±: 

Proposition 18 Let X n be as in If we consider random parameters ao < s t 
S , ai < st Si and Pi < st fix, then \X n \ < st X n , X n < st X n and X n < cx X n . 



Proof. Since a n and cr^are increasing functions of the parameters, if ao < 



si 



ocq, ax <st ol\ and Pi < s t Pi it follows that a n <«t o~ n and a n < s t <J„. As in the 



X n 



and X? <,* X"i. From Lemma 



proof of Theorem [4] it follows that \X n \ < st 

0o"„ <st a n implies that X n ^. cx X n . m 

The last point is to prove the convex comparison of the sums S n ; again, this 
is nontrivial since the X n are not independent; we provide a proof in the case 
of symmetric innovations. 

Theorem 19 Let X n be as in ([?[) and S n as in {5p. Let Si be symmetric. If 
we consider random parameters ao <s* So, ai < s t Si and /3i < s t Pi, then 

Sn — cx S n . 

Proof. As before, we write 

S n = o-q£ Q + gi{s , ...,a ,ai,/3i)ei + ... + g n (e , a , ai,fii)e n 

where the functions gi are nondecreasing in the parameters ao, o>i, pi. Let <f> be 
any convex function. We first want to prove that E [<p (S n )] is nondecreasing 



13 



in the parameters cto,ai,0i. From the symmetry of the innovations £j we can 
write: 

E[<f>(S n )] = S eo! ... !£ „[^(o- O e + ••• + g n e„)} 

= Ee ,...,e n [ ^ ^(o-oPoEO + — +f>nflVi£n)l{ Eo >0 e n >0}]- 

peP„+i 

Denoting by 

h(ao,...,ao,ai,p 1 )= 4>(<?oPo£o + ■■■ + Pn9nS n ), 

pePn+i 

we see that h is nondecreasing in ao,cti,/3i; indeed we can compute: 

cfti 

t, — = <i> (o- p e + ■■• + p n g n £ n ) ■ (pi£ig[ + ■■■+ Pn£ n g'n) > 

peP n +i 

from the multivariate covariance inequality, as in the proof of Lemma 1101 The 
same reasoning shows that > and > 0. 

It follows that E [4> (S n )] is is nondecreasing in ao, a±,0i; but then if ao < s t ot-o, 
ai < st ai and f3\ < st Pi, 

E[4>{S n {a ,ai,P- i )] < E[cj)(S n {a ,ai,Pi)] 

that is S n ~^ C x S n . ■ 
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